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ABSTRACT 

We investigate the quantity and composition of unburned material in the 
outer layers of three normal Type la supernovae (SNe la): 2000dn, 2002cr and 
2004bw. Pristine matter from a white dwarf progenitor is expected to be a 
mixture of oxygen and carbon in approximately equal abundance. Using near- 
infrared (NIR, 0.7 — 2.5 jum) spectra, we find that oxygen is abundant while 
carbon is severely depleted with low upper limits in the outer third of the ejected 
mass. Strong features from the O I line at X rest = 0.7773 fim are observed through 
a wide range of expansion velocities ~ 9 — 18 x 10 3 km s _1 . This large velocity 
domain corresponds to a physical region of the supernova with a large radial 
depth. We show that the ionization of C and O will be substantially the same in 
this region. C I lines in the NIR are expected to be 7 — 50 times stronger than 
those from O I but there is only marginal evidence of C I in the spectra and none 
of C II. We deduce that for these three normal SNe la, oxygen is more abundant 
than carbon by factors of 10 2 — 10 3 . Mg II is also detected in a velocity range 
similar to that of O I. The presence of O and Mg combined with the absence of C 
indicates that for these SNe la, nuclear burning has reached all but the extreme 
outer layers; any unburned material must have expansion velocities greater than 
18 x 10 3 km s _1 . This result favors deflagration to detonation transition (DD) 
models over pure deflagration models for SNe la. 

Subject headings: infrared: stars — line: formation — line: identification — supernovae: 
general 
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1. Introduction 

There is general agreement that Type la supernovae (SNe la) result from the combustion 
of a degenerate carbon/oxygen (C/O) white dwarf (WD) star as first predicted by Hoyle and 
Fowler (1960). The currently accepted scenario for SNe la is the explosion of a C/O WD 
near the Chandrasekhar mass (M C h ~ 1.44M Q ). The explosion is triggered by compressional 
heating near the center of the WD. 

A C/O WD is the final evolutionary stage for stars with a main sequence mass of 
« 3 — 8M . Initially the progenitor has insufficient mass to produce a SN la. To reach 
Mch, the WD must accrete additional matter through Roche-lobe overflow from a main 
sequence, red giant, or WD companion. Quasi-steady shell burning of accreted material in 
conditions of low pressure and high temperature produces a mixture of C and O with a mass 
ratio near unity (Dominguez, Hoflich & Straniero 2001). If two WDs are entirely mixed in 
a merger scenario (Webbink 1984), the C to O ratio by mass will also be close to 1:1 (Benz, 
Thielmann, & Hills 1989). 

One of the keys to understanding the explosion is the propagation of the nuclear flame. 
Two modes of burning can be distinguished: detonation, with burning velocities slightly 
above the speed of sound, and deflagration, with burning front velocities well below sound 
speed. The two models for SNe la currently under discussion are Pure Deflagration (Ivanova, 
Imshennik, & Chechotkin 1974; Nomoto, Sugimoto, & Neo 1976) and Delayed Detonation 
Transition (DD) (Khokhlov 1991). Despite some differences of the numerical treatment of 
nuclear burning fronts by different groups, both models share the same initial sequence: 
nuclear burning begins near the center of the WD and the burning front becomes Rayleigh- 
Taylor unstable. These instabilities result in large scale mixing of burned and unburned 
matter. During the early burning phase, energy is deposited into the WD causing the outer 
layers to expand with velocities close to the speed of sound (Hoflich & Khokhlov 1996). 
Classical deflagration models cease burning in this phase and leave a chemical structure 
characterized by a thick layer of unburned C and O surrounding inner regions of intermediate- 
mass and iron-peak elements. 

Recent refinements with multi-dimensional deflagration models predict that, within the 
context of standard models, a substantial mass of unburned material is mixed toward the 
center of the SN and burning products are mixed outward toward the surface. Primordial and 
processed elements intermingle throughout the radial chemical structure rather than forming 
distinct concentric layers (Gamezo, et al. 2003; Ropke and Hillebrandt 2005). Significant 
quantities of material remain unburned because the deflagration burning front is unable to 
catch up to the expansion of the WD near the surface (Hoflich & Khokhlov 1996; Khokhlov 
2001). Thus the outer layers in current deflagration models remain unburned by means of 
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causality and large pockets of unburned C/0 will be found at various depths because of 
the nature of the Rayleigh- Taylor instabilities. These generic properties of pure deflagration 
models seem to be in contradiction to observations of so called "Branch-normal" SNe la 
(Branch, Fisher, & Nugent 1993) that require vanishing amounts of unburned matter and 
radially layered chemical structures (Barbon et al. 1990; Hbflich 1995; Wheeler et al. 1998; 
Fisher et al. 1998; Hamuy, et al. 2002; Marion et al. 2003). 

DD models also begin with subsonic deflagration burning, but they make a transition to 
a supersonic detonation front. For models of normal SNe la, the detonation eradicates most 
of the chemical inhomogeneities, producing a radially stratified chemical structure with very 
little unburned matter left at the surface. We note that questions remain about how the 
deflagration to detonation transition occurs. For recent reviews of SN la physical models see 
Branch (1999), Hillebrandt & Niemeyer (2000), and Hdflich (2005). 

In this paper, we investigate the presence of unburned material in the outer layers of 
three normal SNe la by examining carbon and oxygen features in near- infrared (NIR, 0.7 
- 2.4 /iin) spectra. The NIR is an excellent region for this search due to the presence of 
several strong lines from C and O ions. We use observational evidence to make estimates 
for the relative abundance of O and C in these SNe la. General methods are used for this 
analysis to produce upper limits on the abundance of these elements and avoid dependence on 
specific models or parameters. Section 2 describes the observations including data acquisition 
and reduction as well as details of the supernovae. Section 3 presents results with line 
identifications and measurements. Section 4 provides methods and calculations for C and O 
locations, ionizations, line strengths, and abundance ratios. Section 5 provides a summary 
and our conclusions. 



2. Observations 

2.1. Data Acquisition and Reduction 

Low and medium resolution NIR spectra were obtained from SNe la using the 3.0 meter 
telescope at the NASA Infrared Telescope Facility (IRTF) with the SpeX medium-resolution 
spectrograph (Rayner et al. 2003). SpeX provides single exposure coverage of the wavelength 
region 0.8 — 2.5 /im. Using a grating and prism cross-dispersers (SXD mode), the spectral 
resolution of the instrument is R — 750 — 2000 and with a single prism (LRS mode) the 
resolution is R = 150 — 250. A serendipitous characteristic of the LRS observing mode is that 
it extends spectral coverage at the blue end of the NIR to about 0.65 //m. The additional 
wavelength region includes O I at \ rest = 0.7773 fim enabling comparison of absorption 
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features from this line to other NIR spectral features. 

Individual exposure times are capped at 150s to limit background noise from OH. A 
typical set is limited to ten exposures for a total of 25 minutes integration time in order to 
remain within timescales of atmospheric variability. The source is nodded along the slit using 
an A_B_B_A_A_B_B_A_A_B pattern. Calibration images are obtained from AOV standard 
stars. Each observation set also includes calibration images from internal lamps for flat 
fielding and wavelength calibrations. 

The data were reduced using a package of IDL routines specifically designed for the 
reduction of SpeX data (Spextool v. 3.2; Cushing, Vacca, & Rayner 2004). These routines 
perform pair subtraction, flat-fielding, aperture definition, spectral tracing and extraction, 
residual sky subtraction, host galaxy subtraction, and wavelength calibration for data ac- 
quired in both the prism mode and the cross-dispersed mode. Corrections for telluric absorp- 
tion were performed using the extracted spectrum of an AOV star and a specially designed 
IDL package developed by Vacca, Cushing & Rayner (2003). These routines generate a 
telluric correction spectrum by comparing the spectrum of an AOV star to a model AOV 
spectrum that has been scaled to the observed magnitude, smoothed to the observed reso- 
lution and shifted to the observed radial velocity. The telluric correction spectrum is then 
shifted to align the telluric absorption features seen in the SN spectrum and divided into the 
target spectrum. 

The individual spectra from each object are combined using using median statistics. 
Noise levels are calculated and recorded by Spextool for each wavelength bin. Because the 
data are subject to bad pixels which can skew the combined mean value and the error on the 
mean, we substitute the median and median absolute deviation (MAD) as robust estimates 
of these statistics. MAD is defined as: 1.4826 x med(\(f\ — med)\) where / = flux, med is 
the median of the combined values, and 1.4826 is a constant introduced on the assumption 
of a Gaussian distribution of initial values. The error on the median at each wavelength is 
equal to the one-sigma noise level and given by M J^ ' , where N is the number of spectra that 
have been combined. 

Detailed procedures to obtain and reduce IRTF data from SNe la are described by 
Marion et al. (2003). 

2.2. The Supernovae 

SNe 2000dn, 2002cr, and 2004bw were identified as Type la from optical spectra, using 
the Si II feature at A rest = 6355 fim which defines the class. The observed characteris- 
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tics of spectral evolution, light curve, and maximum luminosity indicate that each of these 
supernovae are, within uncertainties, Branch-normal SN la events. 

Spectra from these objects were obtained using the LRS mode of the SpeX instrument at 
R = 250. Each of the spectra contain a prominent feature from the O I line at A rest = 0.7773 
/im. The spectrum from SN 2000dn is found in Marion et al. (2003) while the spectra from 
SNe la 2002cr and 2004bw have not previously been published. 

SN 2000dn was observed on October 2, 2000. A relative B-band light curve for SN 
2000dn provided by W. D. Li and M. Papenkova of the Lick Observatory Supernova Search 
team indicates that B max = October 6. We estimate the date of V max to be B max + 2 days, 
making V max = October 8. The recession velocity of the host galaxy (IC 1468) reported 
by LEDA is 9599 km s _1 which makes it a very distant object to observe with SpeX. The 
estimated My = —19.0 mag is slightly dimmer than the canonical value for SNe la of -19.3 
mag, but the decline rate of 0.8 (Marion et al. 2003) suggests that the luminosity estimate 
may be low. This appears to be a normal SN la near the lower end of the decline rate 
distribution. 

SN 2002cr was observed on May 8 & 10, 2002, under adverse weather conditions. The 
spectrum contains data only from May 10 because the May 8 data is not useful below 0.8 
/im. Matheson et al. (2002) estimate that V max = May 8 from optical spectra. A second 
estimate by Meikle, Mack & Marcha (2002) using optical spectra finds the date of V max to be 
May 9 or 10. The Sternburg Supernova Catalog lists V max as May 13 without reference. The 
development of Mg II features in this spectrum is consistent with other NIR spectra obtained 
a few days before V max from normal SNe la with well-documented light curves. NIR spectral 
development favors the May 10 date for V max . The recession velocity of the host galaxy (NGC 
5468) reported by LEDA is 2875 km s -1 . Amateur photometry indicates My ps —19.1 mag 
which is consistent with a normal SN la. We note that this spectrum appears to be shifted 
in wavelength space from other spectra in our sample. We have reviewed our wavelength 
calibrations and conclude that this shift is real and is an example of diversity in normal SNe 
la. 

SN 2004bw was observed on May 31, 2004. The LRS observing mode was chosen for 
this object to maximize signal-to- noise. Foley, et al. (2004) estimate the date of V max to 
be June 6. The NIR spectral development is characteristic of SNE la about a week before 
maximum which agrees with the estimate by Foley, et al.. The recession velocity of the host 
galaxy (MCG +00-38-19) is reported by LEDA to be 6418 km s -1 . Photometry is sparse for 
this object and amateur observations provide a wide range of possible values for My from 
-18.8 to -19.5 mag. Despite the sparse data, all indicators (optical and NIR spectra plus My 
estimates) agree that this object is a normal SN la with a detached high velocity component 
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in the NIR Ca II triplet (see discussion §3.6). 

3. Results 
3.1. The Spectra 

LRS spectra from SNe la 2000dn, 2002cr, and 2004bw are presented in Figure 1. They 
are normalized to unity at 1.0 /im and shifted by a constant for clarity. The Doppler shifted 
positions for O I at 0.7773 /im, the Ca II "infrared triplet" with mean value of 0.8579 jum 
and C I at 1.0693 /im are marked in the figure with vertical dotted lines at 11,000 km s _1 . 
The shifts for Mg II at 0.9227 and 1.0927 /im are marked at 12,000 km s -1 . O I, Ca II, and 
Mg II produce strong absorption features in these spectra, but there is no evidence for C I. 

Superimposed on the raw spectra in Figure 1 are the same spectra after smoothing with 
a Fourier Transform. A detailed description of the procedure and parameters for smoothing 
are found in Marion et al. (2006, in prep.). The smoothed spectra for these three SNe la are 
used for all subsequent figures in this paper. 

Measured Doppler shifts are used to estimate expansion velocities. In regions where the 
lines are formed close to the photosphere, as they are near the absorption minima, projection 
effects make the Doppler shift estimates lower limits on expansion velocities (§4.1). 

We note that the spectrum from SN 2002cr is shifted slightly in wavelength relative 
to the other spectra. This spectrum has Doppler velocities that are 1,000 - 1,500 km s _1 
lower than other SNe la in our sample. In all other respects, this spectrum is similar to the 
others indicating that the chemical structure is approximately the same. This spectrum was 
obtained at an epoch 5-6 days after the two spectra so the Doppler velocities are expected 
to be slightly lower. Small velocity differences are indications of diversity within the class of 
normal SNe la. 

Inspection of Figure 1 shows that the energy flux from SNe la drops off rapidly with 
increasing wavelength. Within a steeply sloping region, the deepest part of an absorption 
feature may not correspond to the wavelength of the maximum departure from the contin- 
uum. This is not a significant effect where the continuum is relatively flat, as in the region 
of the Si II line at 6355 /im, but it can potentially affect features in the NIR by ~ 0.005 /im 
or pa 1,500 km s _1 . To remove this influence, we fit a local continuum in the region of each 
absorption feature and normalize the spectrum to a flat continuum before identifying the 
wavelength of absorption minimum. 

Figure 2 shows the same three LRS spectra along with six SXD spectra from normal 
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SNe la to demonstrate the similarities in spectra from normal SNe la at similar epochs. The 
SXD spectra in this figure have not been smoothed. The spectra are normalized as in Figure 
1 and the wavelength region is reduced to enlarge the region under discussion. The LRS 
spectra are conspicuous due to their high S/N and the fact that they extend to wavelengths 
shorter than 0.8 jum. All nine spectra exhibit strong absorption features from Ca II observed 
near 0.82 /im and Mg II observed near 0.89 and 1.05 /im. The LRS spectra show O I near 
0.75 /im. There is no obvious detection of C I in the spectra. 

Six of the nine spectra in Figure 2 were published earlier by Marion et al. (2003). The 
exceptions are the spectrum from SN 2002cr and both spectra from SN 2005am (obtained 
at -5 and -0 days before V max ) which have not been previously published. 

3.2. O I Features 

The locations of the four strongest O I lines in the NIR are displayed in Figure 3 at a 
Doppler velocity of 11,000 km s -1 . To help with line comparisons, a condensed version of 
Tabel 2 (§4.3) gives rest wavelengths, shifted wavelengths and line strengths normalized to 
the strongest line from each ion. A strong absorption feature from the 0.7773 /im line is 
found near 0.75 /im in each of the LRS spectra. This is the only O I line detected in these 
spectra, but the other lines are expected to be significantly weaker with estimates of 30%, 
6%, and 2% of the 0.7773 /im line strength. In addition, the lines at 0.8446 and 0.9264 /im 
are close to much stronger features from Ca II and Mg II respectively. The non-detection of 
O I lines other than 0.7773 /im is easily explained by weak lines and obscuration by other 
features. 

Atmospheric absorption from O I at 0.7773 /im will contaminate the signal from the 
SNe. Figure 4 shows that the size of the telluric feature is much smaller than the observed 
absorption in the supernova spectrum. Wavelength in this figure is as observed and not 
shifted to the rest wavelength of the host galaxy as in the other figures. The spectra in 
Figure 4 are normalized to 1.0 at 0.745 /xm. The absorption feature in the calibration star 
(dotted line) is clearly smaller than the absorption from the SN (solid line). The telluric 
correction spectrum is the dashed line and the final reduced and smoothed spectrum is the 
thick solid line. The absorption feature in the final spectrum appears to be weaker in this 
figure than it appears in other figures because the wavelength region is small compared to 
the flux region. 

The absorption features from O I at 0.7773 /im are expanded in Figure 5. The top 
panel shows line profiles from the LRS spectra with estimated locations for the continua 
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indicated by dashed lines which are colored to match the spectra. The sizes and shapes 
of these features are generally similar with the most notable difference being the truncated 
bottom in the spectrum of SN 2000dn. The shallowness of this feature is likely due to noise. 
One-sigma noise in the region of this feature is 0.07 — 0.10 for SN 2000dn, 0.04 — 0.08 for 
SN 2002cr, and 0.04 - 0.06 for SN 2004bw. 

In the lower panel, the same line profiles are normalized to a flat continuum and plotted 
in velocity space rather than wavelength space. Expansion velocities are determined by the 
Doppler shift from the rest position of the O I line. The Doppler velocities of the absorption 
minima range from 9 — 11 x 10 3 km s _1 with the blue wings crossing the detection limit at 
18 — 20 x 10 3 km s _1 . Thus the depth of the line forming region in velocity space is ~ 9 x 10 3 
km s -1 for O I in these spectra. 

3.3. C I Features 

The locations the five strongest C I lines in this region are plotted in Figure 3 at a 
Doppler velocity of 11,000 km s _1 . The strongest C I line in the NIR at X re st = 1-0693 
/im is expected to be almost 50 times stronger than the O I line at 0.7773 /im (given equal 
abundance and departure from LTE, Tabel 2), but there are no strong features from this 
line in Figure 3 or in the nine spectra displayed in Figure 2. A few possible weak detections 
seen in Figure 3 are actually smaller than one-sigma noise. The best detection candidate 
is the feature in the spectrum of SN 2004bw found near 1.135 /im and possibly due to the 
1.1756 /im line at an expansion velocity of 10,500 km s _1 . 

Figure 6 shows the spectrum from SN 2004bw in the wavelength range 0.98 — 1.18 /im 
with C I lines from 1.0693 and 1.1756 jum marked at 10,500 km s -1 . We use this figure 
to investigate both the most likely C I detection and the strongest C I line in our data. 
The large absorption feature with a minimum near 1.05 /im is due to the Mg II triplet with 
mean rest wavelength of 1.0927 /im. The spectrum is plotted as a solid line, the estimated 
continuum is plotted as a dashed line, and the one-sigma noise level is plotted with dash- 
dot lines. Noise levels increase in region 1.08 — 1.15 /im due to a reduction in atmospheric 
transmission by about 25%. The actual region of increased opacity is from m 1.1 — 1.16 /im, 
but in the figure the entire spectrum is shifted to the rest wave length of the host galaxy. 

We note that the estimated line strength of the 1.0693 /im line is 7.2 times stronger 
than the line strength of the 1.1756 /im line. If the feature found near 1.135 /im is due to C I 
at 1.1756 /tm producing a line depth of 0.045, then the line from C I at 1.0693 /im should 
produce an absorption feature with a depth of 0.324. A synthetic feature of this size appears 
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in the figure as a dotted line with absorption minimum near 1.03 jum. The presence of such 
a feature is clearly ruled out and we infer that the weak feature at 1.135 /im is not due to 
CI. 

We can extend this argument to all C I lines in Figure 3. If any of the other possible 
detections of other C I lines were real, then the 1.0693 /im line would be expected to produce 
a prominent feature. No such feature is observed in any of the spectra in our sample. In 
addition, the Doppler velocities of the possible detections within individual spectra vary by 
2 — 5 x 10 3 km s _1 . We conclude that none of these features are due to C I. 

Neutral carbon has been detected in subluminous SN la 1999by (Hoflich et al. 2002). 
The presence of carbon in this supernova is interpreted to indicate that burning took place 
at lower temperatures than for normally bright SNe la. Gerady, et al. (2004b) also find C I 
in NIR spectra from Type Ib/c supernovae. These reports indicate that neutral carbon can 
be easily detected in spectra from SNe la. They reinforce the conclusion that the absence of 
C I features in our data is significant. 

3.4. C II and O II 

There is no evidence for C II in these spectra. Branch, et al. (2003) and Garavini, et 
al. (2005) report the identification of C II features from the line at 0.6580 /im in optical 
spectra of normal SNe la. Unfortunately that wavelength region is not accessible with the 
spectra in our sample. There are few C II lines in the NIR and the minimum excitation 
levels are much higher than for C I (Tabel 2). Lines from different ionization stages cannot 
be directly compared, but these very high excitation levels suggest that C II lines will require 
significantly greater abundances than C I to produce detectable features. Even at the upper 
limit of expected temperatures, C II lines in the NIR remain significantly weaker than the 
0.6580 fim line for which the detection is reported. 

There is no evidence in the spectra for O II. Minimum excitation levels for the O II 
lines are greater than 25 eV which is well above the 9 — 11 eV for the O I lines and make it 
unlikely that a detectable signal will be produced in the region covered by these spectra. 

3.5. Mg II Features 

The strongest NIR lines from Mg II are 0.9229 and 1.0927 /im. All spectra in Figure 2 
exhibit absorption features from both lines, near 0.89 and 1.05 /im respectively, with Doppler 
velocities near 12,000 km s _1 . No convincing detections are made for Mg II from lines at 
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0.7890, 0.8228, or 1.0092 fim. 

The top panel of Figure 7 shows the LRS spectra in the region of the absorption feature 
produced by the Mg II line at 1.0927 /im. The estimated location of the continuum for each 
spectrum is marked by a dashed line in the same color as the corresponding spectrum. The 
profiles from this line are generally similar, with the spectrum from SN 2000dn displaying a 
double bottom that is probably due to noise. One-sigma noise in this region is 0.02 — 0.03 
for SN 2000dn, 0.005 - 0.01 for SN 2002cr, and 0.01 - 0.015 for SN 2004bw. 

The lower panel shows the same Mg II features plotted in velocity space and normalized 
to a flat continuum. The velocity minima are in the range 10 - 13 x 10 3 km s -1 and the blue 
wings reach 20 — 21 x 10 3 km s _1 . The depth of the line forming region for Mg II in these 
spectra is ~ 10 x 10 3 km s _1 which is large enough to make nearly all of the magnesium 
singly ionized (§4.2). The velocity range for Mg II is very close to the velocity range we 
identified for O I. 



3.6. Ca II Features 

The spectral region 0.7 — 1.2 fim is dominated by a huge P Cygni feature produced by 
the Ca II infrared triplet with a weighted mean rest wavelength of 0.8579 /zm. The very low 
excitation value of this line (1.7 eV, Table 3) makes it strong enough to obscure the signal 
from most nearby lines. Doppler velocities of the absorption minima are between 11,000 and 
13,000 km s -1 . 

Many SNe la show a detached high velocity component in the absorption feature pro- 
duced by the Ca II IR triplet (Hatano, et al. 1999; Li et al. 2001; Wang et al. 2003; Mazzali 
et al. 2005a,b; Quimby, et al. 2006). The spectrum of SN 2004bw clearly shows a detached 
high velocity Ca II feature at ps 22,000 km s -1 (Figure 3). Gerardy et al. (2004a) have have 
argued that this high velocity feature is the spectral signature of a region where the shock 
from the explosion encounters and compacts circumstellar material forming a dense shell. 
The shell is composed of material from the circumstellar environment, possibly including an 
accretion disk, as well as the extreme outer layers of the supernova ejecta. The disk has a 
high scale height giving it a large covering factor that makes it observable from a wide range 
of angles (Quimby, et al. 2006). Gerardy et al. provide a figure showing the mass of the shell 
as a function of expansion velocity. Our measurement of v exp ~ 22, 000 km s _1 corresponds 
to ps 2 x 1O _2 M for the mass of the shell. 
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4. Limits on Carbon Abundance 

4.1. The Location of Unburned Material 

If the features observed in these spectra are due to unburned material, then C should 
occupy the same physical space. Alternatively, oxygen can be produced by explosive carbon 
burning in regions where the temperature was sufficient to burn C but not hot enough to 
consume O. In that case, magnesium is expected in the same physical region as O because 
Mg is also a product of carbon burning. Since oxygen-rich domains may exist in SNe la from 
both burned and unburned sources, it is possible that C will not be found in all areas where 
O is detected. 

Oxygen has been detected in these spectra with Doppler velocities of ~ 9 — 18 x 10 3 km 
s -1 , measured between the absorption minimum and the blue wing at the detection limit. 
The full range of Doppler velocities is required to identify the entire line forming area and 
not to be restricted to the region closest to the photosphere. The red half of an absorption 
feature does not provide useful information because it is contaminated by emission. All lines 
contain an emission component due to a non-vanishing line source function. 

Expansion velocity in SNe la is directly proportional to radial distance as a result of 
homologous expansion. A continuous range of expansion velocities can therefore be used to 
represent contiguous physical space in the supernova. Velocities uncertainties from spectral 
resolution are ps ±1, 200 km s" 1 but Doppler velocities can be determined to « ±500 km 
s _1 by interpolating between pixels. 

We note that the continuum forming region is extended. If the line-formation region for 
a specific ion is close to this photosphere, many incoming photons are absorbed from matter 
with a projected velocity of v x sin(9) where 9 < n/2. This projection effect shifts the 
observed position of absorption minima to the red which causes the true expansion velocity to 
be underestimated by 10 — 20%. An example is the O I line-forming at the epochs considered 
here. For lines formed far outside in the last scattering radius for continuum photons, the 
photosphere appears to be more point-like and the projection effect is minimized. Typical 
examples are the NIR Mg II lines at later times when the Mg-rich layers are well above the 
photosphere. 

An estimate of the relationship between expansion velocity and mass distribution has 
been calculated for DD models by Hoflich et al. (2002). All DD models release about the 
same amount of energy because nearly the entire C/O WD is burned, so the velocity/mass 
relationship is not sensitive to specific model parameters. Hoflich et al. determine that about 
half of the mass in a SN la is exterior to a velocity of 9,000 km s _1 but only about 5% of 
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the mass is beyond 18,000 km s _1 . This means that the line-forming region for O I in the 
SNe la we have observed extends to the extreme outer layers and includes more than one 
third of the total mass of the supernova. The lack of carbon in this region implies that 5% 
is the upper limit on the mass remaining unburned with C and O in approximately equal 
abunance. 

This estimate, that ~5% of the mass has a velocity exceeding 18,000 km s _1 , is an upper 
limit based on the assumption that nearly all of a C/O WD progenitor has been burned. If 
a different model is used, for instance a pure deflagration model, then less mass is burned 
and less energy is released. The reduction in total kinetic energy implies that even less mass 
will have expansion velocities exceeding 18,000 km s _1 . 

4.2. Finding the Dominant Ionization 

In SNe la, transitions from one ionization stage to another are rapid in both time 
and velocity space due to the dominance of radiative transitions over collisional processes. 
Model calculations for radiative transfer in SNe la envelopes indicate that UV opacities in 
the ejecta are dominated by Fe-group lines that create a large optical depth for UV photons. 
The distribution of the photo-ionization spectrum is determined by Fe opacities that generate 
a smooth quasi- continuum. The ionization balance for other elements is hardly dependent 
on their own photo-ionization boundaries. As a result, elements with similar ionization 
potentials cover similar spatial and velocity ranges in models for SNe la, with boundaries 
defined by Fe transitions (Hoflich 1995). 

Table 1 shows that the first ionization potential for carbon is 11.26030 eV and for 
oxygen, 13.61806 eV. These values are between the first ionization (7.9024 eV) and the 
second ionization (16.1878 eV) of iron. The energy required to doubly ionize both C and O 
is greater than the ionization edge for Fe III. Fe opacities restrict both carbon and oxygen to 
be primarily neutral in the physical zone bounded by Fe ionizations at 7.9 eV and 16.1 eV. 

Small differences in the ionization potential between two elements will create a narrow 
band in radial space where one element is ionized and the other is not. The size of this 
region is model dependent but it is on the order of 1 — 2 x 10 3 km s _1 (Hoflich 1995) which is 
much smaller than the observed depth of the O I line forming region pa 9 x 10 3 km s _1 . This 
large radial depth eliminates the possibility that within the line forming region a significant 
amount of carbon will exist in an ionization state different from oxygen. 

Elements with a first ionization potential similar to magnesium (7.64624 eV) will also 
be singly ionized in the velocity range observed for Mg II. Table 1 shows atoms that fit that 
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description include: Si, Ca, Ti, Mn, Fe, Co, and Ni. 

4.3. Estimating Line Strength in LTE 

Level populations in SNe la are far from local thermodynamic equilibrium (LTE), how- 
ever the line strength of transitions from elements with similar ionization and excitation 
levels show similar departures from LTE so they can be directly compared (Hoflich 1995). 
Uncertainties due to differences in the departure coefficients are a factor of 2-3 which does not 
alter our conclusions. The use of LTE for calculating line strengths removes the dependence 
of our estimates on specific SN la model parameters or methods. 

Temperatures in the line producing region at this epoch are expected to be between 
5,000 - 10,000K (Hoflich, Mueller & Khokhlov 1993; Khokhlov, Mueller & Hdflich 1993). 
Line strength is calculated for each line at temperatures of 5,000K and 10,000K using (gf) * 
lO-eW where 6 = 5040/T. The results are given in Tables 2 and 3. To facilitate comparison 
of line strengths, the computed values at 5,000K are normalized to the strongest line for each 
ion. Rest wavelengths are given in air. Many of the indicated wavelengths are a mean value, 
weighted by oscillator strength, for two or more lines. For these multiplets, the log(g/) value 
and excitation minima are listed for the strongest line in the blend. 

Temperature estimates can be derived from the line strength estimates in combination 
with observations of Mg II features. The spectra in Figure 2 show strong features from 
Mg II lines at 0.9229 and 1.0927 //m but no detections from lines at 0.7890, 0.8228, or 1.0092 
/im. These results are consistent with relative line strengths calculated for a temperature 
of 5,000K. Calculations for the same lines using a temperature of 10,000K suggest that the 
0.7890 /im line will be slightly stronger than the 1.0927 /im line and that the 1.0092 line 
should be detectable with more than 40% of the strength of the 1.0927 /im line. These line 
strengths are not consistent with the spectral features. We conclude that the region probed 
by NIR spectra before maximum light is closer to 5,000K than to 10,000K. 

4.4. Calculating Relative Abundance Ratios for Carbon and Oxygen 

The abundance ratio is the quotient of the line depth ratio divided by the line strength 
ratio. O to C abundance ratios are determined separately for each of he five strongest C I 
lines covered by our observations in each LRS spectrum. (C I at 1.4543 /im is not included 
because it falls in a wavelength region with poor atmospheric transmission). Line depth 
ratios are calculated at all velocities in the range of the O I line-forming region (9 — 18 x 10 3 
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km s" 1 ). The flux array from O I at 0.7773 /im is divided by the flux array for each of the 
C I lines. Wavelength intervals are not constant throughout the spectra, so we interpolate 
between C I velocity data points to match the velocity spacing of the O I data. The line 
strength ratio for each pair of lines is the quotient of the estimated line strength for the O I 
line divided by the estimated line strength for the C I line (Tabel 2, T = 5,000K). Thus we 
determine the O to C abundance ratios over a range of velocities corresponding to a large 
radial depth. 

We illustrate this procedure by calculating the O to C ratio for the three LRS spectra 
using the strongest lines from each ion at a specific velocity. We measure the maximum 
line depths of the O I lines and determine the velocities at maximum depth from the data 
displayed in Figure 5. The depths of the C I 1.0693 /xm lines are measured in each spectrum 
at the same velocities as for O I. The C I line depths are found to be less than the one- 
sigma noise at those locations, so we substitute the noise values for the C I line depths. 
This substitution makes the O/C line depth ratios lower limits. Table 4 gives the measured 
velocities, line depths, line depth ratios, estimated line strengths at 5,000K, the line strength 
ratios, and lower limits for the O/C ratios. Lower limits for the oxygen to carbon abundance 
ratios at the velocities of O I maximum depths are found to be 4,994 for SN 2004bw, 957 for 
SN 2000dn and 2,304 for SN 2002cr. 

Figure 8 shows estimated oxygen to carbon abundance ratios as a function of Doppler 
velocities. O to C ratios have values of 10 2 — 10 3 near the maximum measured depths for O I 
and approach unity near 18,000 km s -1 where both O and C abundances become comparable 
to the one-sigma noise. In many cases, such as the C I line at 1.0693 /im in the spectrum of 
SN 2004bw, the O/C ratio substantially exceeds unity to velocities beyond 20,000 km s _1 . 
At nearly all velocities, the indicated abundance ratios are lower limits because no signal is 
detected from carbon greater than the one-sigma noise level. 

The top panel in Figure 8 displays ratios for O I to C I at 1.0693 /im, the center panel 
gives ratios for C I at 0.9087 and 0.9406 /im, and the lower panel presents ratios for C I at 
0.9639 and 1.1756 /im. Rest wavelengths of the C I lines and their estimated line strengths 
are given in the top right corner of each panel. The slopes plotted in the figure do not always 
follow a monotonically descending path in direct proportion to the O I line depth. Instead, 
the shape of each plot responds to variations in local noise levels and to locations where the 
spectrum in the velocity region of a C I line dips below the one-sigma noise. 
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5. Discussion 

The amount of material remaining unburned after the explosion is an important dis- 
criminant between the leading models for SNe la. Deflagration models predict that large 
quantities of unburned material will remain in the outer layers of the SN la after the ex- 
plosion while Deflagration to Detonation Transition (DD) models predict that essentially 
all progenitor matter will experience nuclear burning. There is agreement that within the 
context of standard models for stellar evolution, the composition of unburned material will 
be carbon and oxygen in approximately equal abundance. Oxygen can be present in SNe la 
from either primordial material or as a product of carbon burning. In unburned matter, C 
and O should both be detectable. In C burning regions, C will have been consumed but Mg 
will be present as well as O. The similarity of ionization potentials for carbon and oxygen, 
the small mean free path for UV photons, and features that cover several thousand km s -1 
imply that carbon and oxygen will be in the same ionization state. 

We have conducted a search for evidence of carbon and oxygen using three NIR spectra 
from normal SNe la. A strong detection is made from O I at 0.7773 /xm with a projected 
velocity range of ~ 9 — 18 x 10 3 km s _1 . The line-forming region comprises about one third 
of the matter in SNe la and O I extends to the outer layers. We do not detect carbon at any 
velocity, in agreement with Barbon et al. (1990); Hoflich (1995); Fisher et al. (1998); Stehle 
et al. (2005). If any matter from these SNe la contains C and O in nearly equal abundance, it 
must have an expansion velocity in excess of 18,000 km s _1 . The upper limit on the mass of 
carbon beyond this velocity is: 3.6 x 1O~ 2 M . We use a DD velocity structure as a limiting 
case for these estimates because alternative models, such as deflagration models, produce 
even less mass above given expansion velocities. 

Strong, broad lines from Mg II are also found in all spectra in our sample. We show 
that the Mg II line-forming region covers a similar velocity range to O I, and hence a similar 
spatial region. This result demonstrates that the domain of explosive carbon burning reaches 
the extreme outer layers of matter probed by these spectra. 

We compare estimated line strengths and measured line depths for O I and C I to derive 
lower limits on the O to C abundance ratios for the SNe la in our sample. We find that the 
oxygen abundance is greater than the carbon abundance by factors of 10 2 — 10 3 at m 11,000 
km s -1 , and remains well above unity to velocities in excess of 18,000 km s _1 . This large 
imbalance and carbon depletion contradict pure deflagration models for standard progenitors 
that predict significant amounts of unburned C and O will remain after the explosion.. 

Patchiness of the carbon and oxygen regions is sometimes invoked as a possible expla- 
nation for the lack of detection of carbon in spectra from SNe la, but the outer layers of 
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published deflagration models are not patchy. The strong O I absorption features in these 
spectra require a large covering factor and cannot be produced by thin clumps of processed 
material. 

Our generalized treatment of the problem has limitations but it adds stability to our 
analysis and eliminates model dependencies. The use of LTE departure coefficients intro- 
duces uncertainty in the estimated line strengths by factors of 2 — 3. The combined uncer- 
tainties in the measurement of expansion velocities are about ±750 km s _1 . The estimated 
location of the local continuum affects measurements of line depths at less than two sigma 
noise levels, but is inconsequential at greater depths. The sum of these uncertainties is not 
nearly large enough to change our results. 

The presence of Mg II in the same region as O I, combined with a lack of detection of 
C I indicate that nuclear burning has taken place in at least the outer one third of the mass 
in these SNe la and extends to more than 18,000 km s -1 . We conclude that for these SNe 
la, carbon is absent at radial depths below a minimum velocity of 18,000 km s _1 because it 
has been explosively burnt to O and Mg. This conclusion is in agreement with DD models 
for SNe la and contradicts pure deflagration models but it does not address the existence of 
C and O in SNe la at velocities below 9,000 km s _1 or above 18,000 km s _1 . 

Quimby, et al. (2006) address the extreme outer layers by showing that Si II is present 
in optical spectra of SNe la to velocities of 24,000 km s -1 . Silicon in this region can 
only be produced from a C/O WD progenitor by explosive oxygen burning at temperatures 
greater than those required to burn carbon into oxygen and magnesium. Non-standard stellar 
evolution processes may be able to enrich oxygen in the outer layers and reduce the amount of 
carbon. This possible explanation for the lack of observed carbon is severely constrained by 
the presence of high temperature burning products in the same region. Carbon abundance in 
the inner regions of the ejecta are constrained by Kozma, et al. (2005) using nebular spectra 
to provide very low limits on the mass of unburned material remaining in the center of SNe 
la after the explosion. 

The results of Quimby, et al. and Kozma, et al. combine nicely with the evidence 
we provide in this paper to demonstrate that nuclear burning must reach every part of 
the progenitor. The lack of carbon detection in all regions of normal SNe la contradicts 
predictions of published 3-D Deflagration models that unburned C will be present at all 
velocities. On the other hand, these observations are perfectly consistent with predictions 
of Deflagration to Detonation Transition models for SNe la. 

Acknowledgments: 

We thank W. D. Li and M. Papenkova for providing light curve data for many of the 



-17- 



SNe in our sample. We thank the individuals at the IRTF for guidance and help with the 
observations. In particular, Alan Tokanaga, John Rayner, Mike Cushing, Bill Golisch, Dave 
Griep, and Paul Sears have been most helpful. We would also like to thank the TAC of the 
IRTF for support and instructive comments. GHM would like to thank Mike Cushing for 
helpful comments. This research is supported in part by NSF grant 0098644 and by NASA 
grant NAG5-7937. CLG is supported through UK PPARC grant PPA/G/S/2003/00040. 

REFERENCES 

Althaus, L. G., Serenelli, A. M., Corsico, A. H., Benvenuto, O. G., 2002 MNRAS, 330, 685 
Barbon R., Benetti S., Cappellaro E., Rosino L., Turatto M. 1990, A&A 237, 79 
Benz, W., Thielemann, F.-K., Hills, J. G. 1989, ApJ, 342, 986 
Branch, D., Fisher, A., and Nugent, P. 1993, A J, 106, 2383 
Branch, D. 1999, ARA&A, 36, 17 

Branch, D., Garnavich, P. Matheson, T., Baron, E., Thomas, R., Hatano, K., Challis, P., 
Jha, S., and Kirshner, R. 2003, ApJ, 126, 1489 

Branch, D. 2004, in: Cosmic Explosions in 3-D, p. 132, Cambridge 

Cusing, M. C, Vacca, W. D. and Rayner, J. T. 2004, PASP 116, 362 

Dommguez, L, Hoflich, P., & Straniero, O. 2001, ApJ557, 126 

Fisher, A., Branch, D., Hoflich, P., Khokhlov, A. 1998, ApJ, 494, 47 

Foley, R. J., Wong, D. S., Morre, M., Filippenko, A. V. 2004, IAU Circ. 8353 

Garavini, G., et al. (The Supernova Cosmology Project) 2005, AJ, in press 

Gamezo, V.N., Khokhlov, A.M., Oran, E.S., Chtchelkanova, A.Y., and Rosenberg, R.O. 
2003, Science, 299, 77 

Gamezo, V.N., Khokhlov, A.M., Oran, E.S. 2005, ApJ, 623, 337 

Gerardy, C. L., HoefTich, P., Quimby, R., Wang, L., Wheeler, J. C, Fesen, R. A., Marion, 
G. H., Nomoto, K, Schaefer, B. E. 2004, ApJ, 607, 391 



Gerardy, C.L., Fesen, R.A., Marion, G.H., Hoflich, P., Wheeler, J.C., Nomoto, K., and 
Motohara, K. 2004, in: Cosmic Explosions in 3-D, p. 57, Cambridge 

Hamuy, M., et al. 2002, AJ, 124, 417 

Hatano, K., Branch, D., Fisher, A., Baron, E., and Filippenko, A. V. 1999, ApJ, 535, 881 

Hoflich P. , Mueller, E., and Khokhlov, A. 1993, A&A, 268, 570 

Hoflich P. 1995, ApJ, 443, 89 

Hoflich, P. and Khokhlov, A. 1996, ApJ, 457, 500 

Hoflich P., Gerardy C, Fesen R., Sakai S. 2002, ApJ, 568, 791 

Hoflich P., Gerardy C, Linder, E., Marion, H. 2003, in: Lecture Notes in Physics, vol. 635, 
p. 203, Springer, Berlin 

Hoflich P. 2005, Nuclear Physics A, in press 

Hillebrandt, W., and Niemeyer, J. 2000 

Hoyle, P., & Fowler, W.A. 1960, ApJ, 1132, 565 

Ivanova, L. N., Imshennik, V. S., Chechotkin, V. M., 1974, Ap&SS, 31, 497 
Khokhlov, A. 1991, ApJ, 245, 114 

Khokhlov, A., Mueller, E., and 1993, Hdflich P. 1993, A&A, 270, 223 

Khokhlov, A. 2001, sgrb.conf. 239, in Supernovae: Models, Explosions, Cambridge University 
Press 

Kozma, C, Fransson, C, Hillebrandt, W., Travaglio, C, Sollerman, J., Reinecke, M., Ropke, 
F. K., and Spyromilio, J. 2005, A&A, 437, 983 

Kurucz, R. L., and Bell, B. 1995, Atomic Line Data, Kurucz CD-ROM No. 23., Cambridge, 
Mass.: Smithsonian Astrophysical Observatory 

Li, W. et al. 2001, PASP 113, 1178 

Marion, G. H., HSflich, P. , Vacca, W. D. and Wheeler, J. C. 2003, ApJ, 591, 316 

Marion, G. H., Hoflich, P. , Gerardy, C. L, Vacca, W. D. Wheeler, J. C, and Robinson, E. 
L., 2006, in preparation 



-19- 

Matheson, T., Jha, S., Challis, P. and Kirshner, R. 2002, IAU Circ. 7891 

Mazzali, P., et al. 2005a, ApJ, 623, 37 

Mazzali, P., et al. 2005b, MNRAS 357, 200 

Meikle, P., Mack, K.-H., and Marcha, M. 2002, IAU Circ. 7891 

Nomoto, K., Sugimoto, D., Neo, S., 1976, Ap&SS, 45, 71 

Quimby, R., Hoflich, P., Kannappan, S. J., Rykoff, E., Rujopakarn, W., Akerlof, C. W., 
Gerardy, C. L. and Wheeler, J. C. 2006, ApJ, 636, 400 

Rayner. J. T., Toomey, D. W., Onaka, P. M., Denault, A. J., Stahlberger, W. E., Vacca, W. 
D., Cushing, M. C. and Wang, S., 2003, PASP 115,362 

Reinecke, M., Hillebrandt, W., and Niemeyer, J. 2002, A&A, 391.1167 

Ropke, F. K. and Hillebrandt, W. 2005, A&A, 431, 635 

Salaris, M., Dominguez, I., Garcia-Berro, E., Hernaz, M., Isern, J., and Mochkovitch, R., 
1997, ApJ, 486, 413 

Stehle, M., Mazzali, P. A., Benetti, S., and Hillebrandt, W. 2005 MNRA, 360, 1231 

Vacca, W. D., Cushing, M. C. and Rayner, J. T. 2003, PASP, 115, 389 

Wang, L., Baade, D., Hoflich, P., Khokhlov, A., Wheeler, J. C, Kasen, D., Nugent, P., 
Perlmutter, S., Fransson, C, Lundqvist, P. 2003, ApJ, 591, 1110 

Webbink, R. F. 1984, ApJ, 277, 355 

Wheeler, J. C, Hoflich, P., Harkness, R. P., Spyromilio, J. 1998, ApJ, 496, 908 



This preprint was prepared with the AAS IATgX macros v5.2. 



-20- 



Table 1. Ionization Potential for SN la Explosion Products 



Atom First (eV) Second (eV) 



c 


11.26030 


24 


,38332 





13.61806 


35 


41730 


Ne 


21.56454 


40. 


.96328 


Mg 


7.64624 


15. 


,03528 


Si 


8.15169 


16 


.34585 


S 


10.36001 


23 


.3379 


Ar 


15.75962 


27. 


.62967 


Ca 


6.11316 


11 


.87172 


Ti 


6.8282 


13 


.5755 


Mn 


7.43402 


15. 


.63999 


Fe 


7.9024 


16 


.1878 


Co 


7.8810 


17. 


.083 


Ni 


7.6398 


18. 


.16884 
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Table 2. Strong Lines for Carbon and Oxygen in the Region 0.6 — 2.4 /zm 



Ion 


A rest (air) 


No. lines 




Excitation 


Line Strength 


Line Strength 


Line Strength 




jum 


in blend a 




Min. (eV) 


5,000K b 


normalized 


10,000k b 


C I 


1.0693 


7 


0.348 


7.48830 


6.307e-08 


1.000 


3.749e-04 


C I 


0.9087 


6 


0.142 


7.48830 


3.925e-08 


0.622 


2.333e-04 


C I 


0.9406 


1 


0.225 


7.68529 


3.008e-08 


0.477 


2.247e-04 


C I 


0.9639 


3 


-0.264 


7.48830 


1.541e-08 


0.244 


9.160e-05 


C I 


1.4543 


1 


-0.110 


7.68529 


1.391e-08 


0.221 


1.039e-04 


C I 


1.1756 


7 


0.661 


8.64774 


8.792e-09 


0.139 


2.007e-04 


C I 


0.8333 


2 


-0.420 


7.68529 


6.811e-09 


0.108 


5.089e-05 


I 


0.7773 


3 


0.324 


9.14671 


1.271e-09 


1.000 


5.177e-05 


I 


0.8446 


3 


0.170 


9.52201 


3.731e-10 


0.294 


2.349e-05 


I 


0.9264 


9 


0.690 


10.74166 


7.285e-ll 


0.057 


1.889e-05 


I 


1.1290 


9 


0.500 


10.98960 


2.645e-ll 


0.021 


9.146e-06 


01 


0.7990 


6 


0.280 


10.98960 


1.594e-ll 


0.013 


5.511e-06 


CII 


0.6580 


2 


0.118 


14.44980 


3.569e-15 


88,393 


6.844e-08 


CII 


0.7235 


3 


0.330 


16.33422 


7.330e-17 


1,815 


1.252e-08 


CII 


1.8905 


2 


0.258 


19.49586 


4.038e-20 


1.000 


2.705e-10 


CII 


1.7846 


3 


0.550 


20.15184 


1.726e-20 


0.427 


2.474e-10 


CII 


0.9903 


3 


1.010 


20.95206 


7.768e-21 


0.192 


2.819e-10 


on 


0.6695 


2 


-0.594 


23.44307 


5.962e-25 


2,155 


3.897e-13 


on 


2.1085 


2 


-1.690 


25.66290 


2.766e-28 


1.000 


2.376e-15 


on 


1.3811 


3 


-1.730 


25.66290 


2.522e-28 


0.912 


2.167e-15 


on 


0.6898 


12 


0.230 


28.70803 


1.960e-29 


0.071 


5.770e-15 



Note. — 1995 Atomic Line Data (R.L. Kurucz and B. Bell) Kurucz CD-ROM No. 23. Cambridge, 
Mass.: Smithsonian Astrophysical Observatory 

a Where the number of lines in the blend exceeds 1, the indicated wavelength is a mean value, 
weighted by oscillator strength, for two or more lines. 

b The estimated line strength is computed using (gf) * 10~ e<ey ) where O = 5040/T. 
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c The normalized line strength is the 5,000K value divided by the strongest line for each species. 



-23- 



Table 3. Strong Lines ofr Mg II and Ca II in the Region 0.6 — 2.4 jj,m 



Ion 


A resi (air) 


No. lines 


log(#/) 


Excitation 


Line Strength 


Line Strength 


Line Strengtl 




(Ltm 


in blend a 




Min. (eV) 


5,000K b 


normalized c 


10,000k b 


Mg II 


0.9227 


2 


0.270 


8.65529 


3.511e-09 


1.000 


8.086e-05 


Mg II 


1.0927 


3 


0.020 


8.86425 


1.216e-09 


0.346 


3.568e-05 


Mg II 


0.7890 


3 


0.650 


10.00000 


3.715e-10 


0.106 


4.074e-05 


Mg II 


0.8228 


2 


0.030 


10.00000 


8.912e-ll 


0.025 


9.772e-06 


Mg II 


1.0092 


3 


1.020 


11.63047 


1.979e-ll 


0.006 


1.440e-05 


Ca II 


0.8579 


3 


-0.362 


1.70005 


8.401e-03 


13,944 


6.042e-02 


Ca II 


1.1876 


2 


0.300 


6.46831 


6.025e-07 


1.000 


1.096e-03 


Ca II 


0.8921 


3 


0.729 


7.05003 


4.193e-07 


0.696 


1.499e-03 


Ca II 


0.8235 


3 


0.621 


7.51535 


l.llle-07 


0.184 


6.812e-04 


Ca II 


0.9906 


2 


0.072 


7.51535 


3.137e-08 


0.052 


1.924e-04 



Note. — 1995 Atomic Line Data (R.L. Kurucz and B. Bell) Kurucz CD-ROM No. 23. Cambridge, 
Mass.: Smithsonian Astrophysical Observatory 

a Where the number of lines in the blend exceeds 1, the indicated wavelength is a mean value, 
weighted by oscillator strength, for two or more lines. 

b The estimated line strength is computed using (gf) * 10~ e<ey ) where O = 5040/T. 

c The normalized line strength is the 5,000K value divided by one of the strongest lines for each 
species. 
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Table 4. Lower Limits on Abundance Ratios for O I at 0.7773 /zm to C I at 1.0693 /iin 



SN 2004bw SN 2000dn SN 2002cr 



Velocity at I max. depth 


10,501 


11,331 


8,447 


I max depth 


0.711 ±0.05 


0.543 ± 0.08 


0.747 ±0.06 


C I depth at same vel. a 


0.007 


0.028 


0.016 


Line depth ratio (^) 


100.5 


19.3 


46.4 


C I line strength estimate 


6.307 x 10~ 8 


6.307 x 10~ 8 


6.307 x 10~ 8 


I line strength estimate 


1.271 x 10~ 9 


1.271 x 10~ 9 


1.271 x 10- 9 


Line strength ratio (^) 


0.02015 


0.02015 


0.02015 


to C abundance ratio b 


4,994 


957 


2,304 



a Measured C I line depth at this velocity does not exceed one-sigma 
noise, so noise values are substituted for line depths. 

b Abundance Ratio = Line Depth Ratio/Line Strength Ratio. 
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Rest Wavelength of Host Galaxy (ju.m) 



Fig. 1. — Three NIR spectra from normal SNe la acquired at the IRTF using the SpeX 
instrument in LRS (single prism) mode, R = 250. The spectra have been normalized to 
unity at 1.0 /mi and shifted by a constant. Superimposed on the raw spectra are the same 
spectra after smoothing with a Fourier Transform. The locations of prominent absorption 
features are identified by a vertical dotted lines at expansion velocities of 12,000 km s -1 for 
the Mg II lines and 11,000 km s _1 for C I, O I, and Ca II. No features are evident from the 
C I line even though it is expected to be ~ 50 times stronger than the O I line. 
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Fig. 2. — Plot of the log of the flux for nine spectra from normal SNe la obtained at or before 
V max . The spectra have been normalized to unity at 1.0 /iin and the vertical displacement 
is arbitrary. The LRS spectra are conspicuous compared to the SXD spectra due to their 
high S/N and the fact that they extend beyond 0.8 /im. The SXD spectra have not been 
smoothed. The locations of prominent absorption features are identified by a vertical dotted 
lines. O I is detected in each of the LRS spectra near 0.75 /im. All spectra in this sample 
contain the Ca II feature observed near 0.82 /im, and the two Mg II features observed near 
0.89 and 1.05 /im. No features are evident from the C I line in any of the spectra even though 
it is expected to be ~ 50 times stronger than the O I line. 
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Fig. 3. — The strongest NIR lines from O I and C I are shifted by 11,000 km s _1 and plotted 
with the LRS spectra. 9,000 km s -1 is ~ 0.01 /mi on the red side of the marked position of 
each line and 18,000 km s -1 is ~ 0.03 /mi on the blue side. To help with line comparisons, 
a condensed version of Tabel 2 is provided in the top right corner. The rest wavelengths, 
Doppler shifted wavelengths at 11,000 km s" 1 , and line strengths normalized for each ion 
are given for strong O I and C I lines. O I from A rest = 0.7773 /mi is found in all spectra 
near 0.75 /im. No reliable detections are made in these spectra for any other lines from O I 
or CI. 
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0.74 0.76 0.78 0.80 

SN 2004bw - Observed Wavelength (fim) 





Fig. 4. — Removal of telluric absorption in the feature produced by O I at 0.7773 fim. All 
spectra are normalized to 1.0 at 0.745 /im. Wavelength in this figure is as observed and not 
shifted to the rest wavelength of the host galaxy as in the other figures. The absorption 
feature in the calibration star (dotted line) is clearly smaller than the absorption from the 
SN (solid line). The telluric correction spectrum is the dashed line and the final reduced 
and smoothed spectrum is the thick solid line. The absorption feature in the final spectrum 
appears to be weaker in this figure than in other figure because the wavelength region is 
small compared to the flux region. 
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Fig. 5. — Enlarged view of the region around the absorption feature from O I at 0.7773 /zm. 
The top panel shows the O I features in three LRS spectra with estimated locations for the 
continua indicated by dashed lines that are colored to match the spectra. In the lower panel, 
the same line profiles are normalized to a flat continuum and are plotted in velocity space 
determined by the Doppler shift from the rest position of the line. The line forming region 
for O I lies between the absorption minimum and the farthest extent of the blue wing. The 
location of the line forming region for O I is similar to that of Mg II (Figure 7). 
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Fig. 6. — The spectrum from SN 2004bw with the location of C I lines from 1.0693 and 
1.1756 /im marked at a velocity of 10,500 km s _1 . The large absorption feature with a 
minimum near 1.05 /im is due to Mg II. The spectrum is plotted as a solid line and the 
estimated location of the continuum is plotted with a dashed line. If the feature found near 
1.135 /im were due to C I, then there should also be an enormous absorption feature near 
1.03 ^m as indicated by the dotted line in the figure. 
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Fig. 7. — Enlarged view of the region around the absorption feature from Mg II at 1.0927 
fim. The top panel shows the Mg II features in three LRS spectra with estimated locations 
for the continua indicated by dashed lines that are colored to match the spectra. In the lower 
panel, the same line profiles are normalized to a flat continuum and are plotted in velocity 
space determined by the Doppler shift from the rest position of the line. The line-forming 
region for Mg II lies between the absorption minimum and the farthest extent of the blue 
wing. The location of the line- forming region for Mg II is similar to that of O I (Figure 5). 
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Fig. 8. — Lower limits for estimated Oxygen to Carbon abundance ratios as a function 
of Doppler velocity. The O I line at 0.7773 /iin is compared to upper limits on the five 
strongest C I lines in the NIR. The top panel displays ratios for O I to C I at 1.0693 /im, the 
center panel at 0.9087 and 0.9406 /im and the lower panel at 0.9639 and 1.1756 /im. Rest 
wavelengths of the C I lines and their estimated line strengths are given in the top right 
corner of each panel. We find that oxygen abundance is greater than carbon abundance by 
factors of 10 2 — 10 3 at ~ 11,000 km s -1 . The lower limit for the O to C ratio remains well 
above unity to velocities in excess of 18,000 km s _1 . 



